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I. INTRODUCTION 



Recently a new approach for calculation the partition function of the 2D Ising 
model on the regular lattice has been suggested |l| . The idea of this approach is as 
follows. An independent generator of Clifford algebra in the matrix representation 
r^, - is assigned to each vertex of the lattice. Thus the number yi of generators of the 
Clifford algebra can not be less than the number of vertexes of the lattice or the 
number of degrees of freedom of the partition function. Some polynomial in these 
generators depending on the "temperature" parameter is defined. This polynomial 
is the matrix of orthogonal rotation in spinor representation in the ^-dimensional 
Euclidean space. The trace of the polynomial is proportional to partition function 
of the 2D Ising model. This statement follows from the fact that the considering 
trace is represented as a sum over all closed self-intersecting |8| contours (loops) on 
the regular planar lattice, and the positive weight depending on the temperature is 
assigned to each edge of the loops. 



In this paper the outlined approach is extended to three-dimensional statistical 
system on on the regular cubic lattice. The sum over self-intersecting surfaces is 
calculated. The closed surfaces are included as well as the surfaces with boundary. 
Each face of the surface is assigned the positive factor it 2 , each edge of the boundary 
is assigned the factor A/i. The real numbers A and fi depend on the temperature. 
Thus, the total positive weight corresponds to each closed surface. However, it turns 
out that the total weights corresponding to the surfaces with boundary depend on 
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the boundary configurations, but not on the surface configurations with the fixed 
boundary. 

The fact that it is impossible to retain only the closed surfaces and to remove the 
surfaces with boundary is the important characteristics of the considered statistical 
sum. This means that the considered system does not involve the 3D Ising model. 
Indeed, the Ising model on three-dimensional regular lattice is dual to the thre- 
dimensional lattice gauge model with the gauge group [2] . The partition function 
of the latter is expressed as the sum over closed self-intersecting surfaces, and each 
face of the surface is assigned the positive factor depending on the temperature 
at that. So it is naturally to interpret the boundary of the surface as creation, 
propagation and annihilation of the fermion pair coupled with the gauge field. It 
will be shown in the subsection 2.5 that these fermions are the lattice analog of the 
Majorana fermions. 

Thus, the calculated partition function can be interpreted as the partition function 
of the gauge system on the regular three-dimensional lattice with the gauge group 
interecting with Majorana fermions |9|]. 



II. FORMULATION OF THE PROBLEM 
A. The Dirac- Clifford algebra in the matrix representation 



The generators of the Clifford algebra {£ x } satisfy the following relations 

6c£y + fw6c = 2 $r,yj x, y,... = 1, ...,/ = 2K, (2.1) 

and they are assumed as hermitian matrixes of the dimension 2 K x 2 . Such di- 
mension is minimally admissible at given number of generators. The theory of such 
matrixes in convenient for physicist form can be found in |3| for The concrete ma- 
trix representation of the Clifford algebra generators is no object. Only some their 
algebraic properties are needed. It follows from the algebra (12.11 ) that the trace of 
any product of the odd number of ^-matrix is equal to zero, and 

£,x£,y = 2 $x,yi tr £,x£,y£,z£,v = 2 ^x,y^z,v ^x,z^y,v ^x,v^y,z) j ( 2-2 ) 

and so on. It is evident that by means of permutations and usage of the equalities 
£| = 1 any product of ^-matrixes is reduced either to the number ±1 (versus the 
number of permutations) or to the product of pairwise different ^-matrixes. Ac- 
cording to Eq. (12.21 ) in the last case the trace of the product is equal to zero, and 
in the firs case it is equal to ±2 . 



3 



B. Description of the degrees of freedom of the system 

Let's enumerate the axes the simple three-dimensional cubic lattice by the numbers 
% = 1, 2, 3. Each vertex t) x has the index composed of three natural numbers 
x = (to, n, /), and 

m = l,...,M, n=l,...,N, I = 1, . . . , L. (2.3) 

For simplicity the numbers M, N, L are considered even. It is assumed that the 
numbers (to, n, I) increase uniformly along the first, the second and the third axis, 
correspondingly. Define the base vectors of the lattice: ei = (1, 0, 0), e2 = (0, 1, 0), 
e3 = (0, 0, 1). Thus, the vector is directed onto positive direction of the z'-th axis. 
Further I believe I = 3MNL. 

Divide the totality of the Clifford algebra generators into 3 groups, each of which 
consisting of MNL generators. Let's designate the first, the second and the third 
group of the generators as {a x }, {/3 X } and {7 X }, correspondingly. The edge con- 
necting the adjacent vertexes o x and t) x +e, is designated as [ x ^ = [ m , n ,z ; j. The face 
with vertexes (x, x + ej, x + e ? - + e^, x + e&) is designated as f x i? % = 1, 2, 3 (here 
j k and the scalar products ej • = • e« = 0). 

By definition, the matrixes a x , /3 X and 7 X are related to the edges [ Xj i, [ x ^ and 
tx,3, correspondingly. Below also the notation 

£x = # = (q x1 7x) , 2 = 1,2,3. (2.4) 

is used. 

C. The definition of the partition function 

Relate to each face the unitary matrix of rotation in the spinor representation: 
The face f x 1 is related with matrix 

U^ nJ) = (A + /X &7x + e 2 ) (A + ^x+e 3 7x) , (2.5) 

the face f 2 is related with matrix 

U (m,n,l) = ( A + ^«x7x+ ei ) (A + ^«x+e 3 7x) , (2.6) 

and the face f x 3 is related with matrix 

U[^ nj) = (A + fia x (3 x+ei ) (A + fia x+e2 (3 x ) . (2.7) 

Here 

ip ib 

x = (to, n, /), A = cos — , /i = sm—. (2.8) 

Zu Zj 
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x+e 2 +e 3 



x+e 3 




x+e ^ s x+ei+e 2 



FIG. 1 



In fig. 1 a fragment of the lattice containing the vertex x and three adjacent faces 
fx, l? fx, 2? fx, 3 is represented, and also the ^-matrixes related to the correspond- 
ing edges are designated. The arrows, curved orthogonally, lying in the faces and 
showing the crossing from one edge to the neighboring edge, represent the rota- 
tion matrixes corresponding to one of the round brackets in the right hand side of 
Eqs. ( 12.51 ). ( 12.61 ) and ( 12.71 ); the directions of the arrows indicate the ordering of 
^-matrixes. For example, the arrows (1) and (2) represent the first and the second 
round brackets in Eq. ( 12.51 ). correspondingly, etc. Since the round brackets in the 
right hand sides of Eqs. ( 12.51 ). ( 12.61 ) and ( 12.71 ) do not commute, so their order is 
significant. Below we need the representation of the quantities ( I2.5I )-( |2T71 ) in the 
form of the following sums 



(m, n, I) 



(A 2 + /i 2 /3 x 7x+e 2 /?x+e 3 7x) + tXfi (/?x7x+e 2 + Ac +e 3 7: 



e 3 ixj 5 



(2.9) 



U (m,n,l) = + ^ 2 «x7x+e i ax+e 3 7x) + tXfJb (a x 7 x+ei + a x +e 3 7x) , 



;2.io) 



u 



(3) 

(m, n, I) 



In the considered theory 



t = 1. 



;2.i2) 



In the case t ^ 1 the matrixes ( I2.9I )-( |2~TT1 ) are not the rotation matrixes and therefore 
the used here calculation method turns out useless. 
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The partition function of the system is proportional to the following trace of the 
ordered product of rotation matrixes in spinor representation 

Z = CtrU, (2.13) 



U = ( P rf) W (0> ( 2 - 14 ) 



M(n,l) — I ^Yl. ) U(m,n,l) ~ ■■■ ^(m-l,n 1 0^(m,n,J)^(m+l,n,0 (2.16) 



U(m,n,l) ^(m, n, l)^(m, n, l)^(m, n, Z)' (^•^'^) 

Here C is some numerical constant which is not interesting here. The symbol 
(P H m ) in ( 12.161 ) denotes the ordered product of the matrices ( 12.171 ) over all mean- 
ing of m at fixed values of the numbers (n, /), so that in this product the matrix 
Uc m ,n,i) is placed to the left of the matrix Ui m > in ,i) f° r m < m'. The ordered products 
in ( 12.151 ) (at fixed (/)) and in ( 12.141 ) are determined analogously. 



D. The closed surfaces with self-intersections 



Let's consider the contribution of zeroth power in the parameter t to the matrix 
( 12.141 ). According to Eqs. ( I2.9I )-( I2TTT1 ) in this case each face gives two contributions 
to U: the first contribution is proportional to the number A 2 and the second one 
is proportional to the matrix {^/, 2 x(the ordered product of four ^-matrixes corre- 
sponding to four edges of the face)}. It is obvious that under the trace operator in 
( 12.131 ) some summands of the matrix IA are inessential. For example, according to 

(E2D 

trC {(A 2 )( 3 ™-V/?x7x + eA+e 3 7x} = 0. (2.18) 

Only that summands of the matrix IA give the nonzero contribution to the partition 
function in which each contained ^-matrix is in even power. It follows from here 
that the considered part of partition function is represented as the sum over closed 
surfaces on the lattice with self-intersections. Each edge of the surface can belong to 
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two or four faces of the surface. In the last case the intersection or self-intersection 
of the surface occur in the edge. To each face f x i of such surface the summand 

from n iy i = 1, 2, 3, proportional to fi 2 and the fourth power of ^-matrix (see 
(12.91 )- ( 12TTT| )). is assigned. The order of the ^-matrix arrangement in the summand 
in U corresponding to the closed surface (not obvious connected) is determined by 
the formulae (I2.13l) -( l2.17l) . It follows from the aforesaid that in the summand in U 
corresponding to any closed surface each ^-matrix can be contained either in the 
power (the corresponding edge does not belong to the surface) either in the power 
2 (the corresponding edge belong to two faces of the surface only) or in the power 
4 (the corresponding edge belong to four faces of the surface). If S is the number 
of faces of a the closed surface, then the contribution to the partition function of 
the surface is equal to (at macroscopic dimensions of the lattice the total number 
of faces of the lattice can be considered equal to 3MNL) 



AZ = C 



2 3 ™/ 2 (A 2 ) 3 ™! (fi/\) 2S , C > 0. (2.19) 



Up to the sign the equality (12.191) is evident. It remains to prove only that each 
closed surface gives the positive contribution to the partition function. 

Consider the contribution to U from the simplest closed surface without self- 
intersections: the elementary cube with the faces f x l3 f x 2 , f x 3, f x + ei i> f x +e 2 ,2? 
fx+e 3 ,3- Here the order of the face enumeration corresponds to the order of con- 
struction from elementary "bricks" the corresponding summand in IA. Thus 

AU = (/?x7x+e 2 /?x+e 3 7x) (ax7x+ ei ax+e 3 7x) (axfeejflx+e^x) X 
x (/5x+ei7x+ei+e 2 /^x+ei+e 3 7x+ei) («x+e 2 7x+ei+e 2 «x+e 2 +e 3 7x+e 2 ) x 

x (,Cyx+e 3 +ei+e 3 Qoc+e 2 +e 3 /5 x +e 3 ) — 1- (2.20) 

In this product each ^-matrix is contained twice. 
To prove the positivity of the expression ( 12.191) in general case let's cut any closed 



surface by the plane parallel to the plane (12), intersecting the middles of the edges 
(m,n,Z;3 with fixed / and m = 1, ... , M, n = 1, . . . , N. Denote this plane by 
In each plane ^ there is the natural structure of the square plane lattice with the 
vertexes designated as t> m ,n which belong to the middles of the edges [ m , n ,Z;3- The 
edges of this lattice are designated by [ TOjn; i and [ m , n; 2- The edge [ m , n ;i connects 
the vertexes t> m ,„ and \) m +i,m an d the edge [ TOj7i; 2 connects the vertexes t> TOj „ and 
b m ,n+i- Thus the edge [ TO>n; i belong to the face f x 2 and the edge [ m ^ n ;2 belong to 
the face f xl . The outlined lattice also is denoted as It is convenient to consider 
that the matrixes 7 TO; „ ; ; with fixed / and m = 1, ... , M, n = 1, . . . , N are related 
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to the vertexes TO;n of the lattice 

Let's fix / and represent by thick segments that edges of the lattice ^ which belong 
to the faces of the considered closed surface. It is evident that the totality of thick 
segments or edges on the lattice ^ forms a closed contour with self-intersections on 
the cubic plane lattice, so that in one vertex o m;n 0, 2 or 4 thick edges can converge. 
It is known that the sum over such closed contours with a certain positive weight 
assigned to each edge is proportional to the partition function of 2D Ising model 
(0], 0])- Let us transform the matrix factors corresponding to the faces f x x and 
fx, 2 (belonging to the closed surface and intersecting with the plane ^) as follows: 

(/?x7x+e 2 /?x+e 3 7x) = (7x7x+e 2 ) (/?x/?x+e 3 ) , ( 2 - 21 ) 

(ax7x+e 1 ax+e 3 7x) = (7x7x+ ei ) (QxQx+e 3 ) • (2.22) 

Note that the matrixes (7x7x+eJ, 2 = 1, 2 commute with all matrixes a x /, (3 X > 
without restrictions and also with the matrixes (7x'7x'+eJ? J = 1, 2 for / ^ I'. 
But the matrixes (7x7x+eJ, i = 1, 2 and (7 x '7x / +e J ), j = 1, 2, for / = generally 
speaking, do not commute since in the set of four matrixes (t*, 7 x +e i? 7x', 7x'+eJ 
can be found coinciding matrixes. It is seen from here and from Eqs. (12.211 ). ( 12.221 ) 



that in the summand in IA corresponding to a closed surface one can remove to 
the left all matrix factors (7x7x+e 2 ) an d (7x7x+ei) without changing the relative 
arrangement of the factors themselves. As a result the sign of the summand does 
not change. Moreover, the product of all removed to the left 7-matrixes is converted 
into 1. Indeed, the ordering product of the removed to the left 7-matrixes (at fixed 
/) corresponds to the ordering product of 7-matrixes in the case of 2D Ising model 
I If . It was proved that in the last case the considered product of 7-matrixes is 
equal to 1. 

Analogously the 7-matrixes for all others / are mutually cancelled. Thus in any 
summand in the matrix IA we can cut out all 7-matrixes without changing the sign. 

Further we work analogously. Let's denote the plane cutting a closed surface 
through the middles of the edges [ mjn ,i;2, m = 1, . . . , M, I = 1, . . . , L by On 
the plane ^ n the natural structure of the right square plane lattice is present. The 
vertexes of this lattice (at fixed n) are designated as tt m ,i- The matrixes Pm^j are 
related to the vertexes \) m j. Taking into account Eqs. ( 12.211 ) and 



(ax/?x+ ei a x +e 2 /?x) = (AcAc+ei) ( a x«x+e 2 ) , (2.23) 

we can (at fixed n) remove from the rest of the ^-matrix product (corresponding to 
the closed surface) all factors (/5 x/ 5 x+ e 3 ) and (/3 x /3 x+ei ) to the left without changing 
the relative arrangement of the factors themselves. As above, the total sign of the 
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^-matrix product does not change, and the product of /^-matrixes at each value 
of n is converted to 1. Further only the product of a-matrixes remains which is 
converted to 1 at each fixed m. 

Thus it is proved that the product of ^-matrixes in each matrix summand in IA 
corresponding to any closed surface is converted to 1, and the formula ( 12.191 ) is true. 



E. The surfaces with boundary 



Let us consider the contributions to the partition function ( 12.131 ) which are nonzero 
relative to the parameter t in ( 12.91 )- ( 12TTTT ) . 

The elementary example of the surfaces with boundary is represented in fig. 2. The 
totality of curved orthogonally thick arrows jjlOj forms the boundary of the surface. 
The surface itself is the shaded face in fig. 2 with the vertexes O2, O3, O4, 06- 

t>5 




FIG. 2 



The contribution to the partition function of the surface with boundary repre- 
sented in fig. 2 is proportional to t 4 . According to the rule the contribution of this 
surface to U is proportional to the expression 

AiW ~ (/?i7i)(«i^i)(o;i72a27i)(^272)(«2^2) = -1. (2.24) 

Here and below the proportionality coefficients are positive. In ( 12.241 ) the factors 
(/3i7i) contained in parentheses are the contributions from the curved arrow in the 
face with vertexes t)i, O2, 06, O7. Generally, the set of factors contained in parenthe- 
ses is the factor in AU either from one arrow in some face (only two ^-matrixes), or 
the factor from the whole face (four ^-matrixes). Each face gives no more than one 
factor which can be either zeroth power relative to the ^-matrixes and proportional 
to A 2 , either the second power (represented by one curved arrow and proportional 
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to t), or the fourth power relative to ^-matrixes. For example, the shaded face 
with vertexes 02, t>3, O4, &6 gives the contribution of the fourth power relative to 
^-matrixes: (a^a^Ti) • 

Now let's consider the contribution to IA from the surface with boundary repre- 
sented in fig. 3. As in the previous example, the curved orthogonally thick arrows 
form the boundary of the surface, and the surface itself consists of the three shaded 
faces. The contribution of this surface to the matrix IA is proportional to t 6 and the 
matrix 

A 2 U ~ (/?i7i)(«i/?i)(«i72a27i)(/?372)(/?273)x 

x(a 2 /3 3 «3/?2)(a374a473)(/?474)(«4/?4) = 1- (2.25) 




FIG. 3 

It is seen from the comparison of the right hand sides of the relations ( 12.241 ) and 
( 12.251 ) that the factor in the summand of the partition function corresponding to 
a surface with boundary can be of any sign. The more detailed consideration of 
this problem shows that the signs of the factors of the surfaces with boundary 
depend on the boundary configuration but not on the surface configuration with 
given boundary. 
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The outlined situation resembles the situation in the case of Majorana fermions on 
three-dimensional lattice coupled with Abelian gauge field. It is shown in Appendix 
B that in the last case the fermion contribution to the partition function also is 
represented as a sum over the surfaces with boundary, and the boundary of the 
surface is interpreted as creation, propagation and annihilation of the fermion pair. 
As in the studied here model, each such surface gives some factor in one of summand 
of all functional integral. The sign of this factor can be either positive or negative, 
and the sign is determined by the boundary configuration but not by the surface 
configuration with given boundary. The modulo of this factor is proportional to e s , 
where e is some positive number and s is the number of faces of the corresponding 
surface. 

Stated above founds to accept the hypotheses that the studied here partition 
function is the partition function of the abelian gauge system with the gauge group 
Z<i (this is evident) coupled with Majorana fermions. Just owing to this hypotheses 
the word " supersymmetric" is introduced into the title of the paper. 

In connection with this hypotheses it is necessary to make the following comment. 

As it is known [5| , it is impossible to construct on the periodic lattice a satisfactory 
variant of the local action of Dirac fermions which transforms into the well known 
Dirac action in the continual limit. In particular, the problem consists in the impos- 
sibility of construction of fermion action on the periodic lattices describing only one 
chiral (Weyl) fermion field and transforming into the usual action of chiral field in 
the continual limit. This problem is known as a "fermion doubling" problem. It is 
known also that at the cost of increase the lattice dimension one can construct only 
one Weyl field on the lattice hypersurface of the codimension one. For example, 
on the 4+1-dimensional lattice one can construct the action of only one Weyl field 
on 4-dimensional sublattice [6]. But this construction does not seem refined and 
therefore prospective. Besides the problem of coupling of the such Weyl fields with 
the gauge field is not solved in general case beyond the framework of perturbation 
theory. 

From here one can make the following conclusion: the approach to the study of 
the lattice gauge systems including fermion degrees of freedom which is based on 
the usage of a lattice action is incorrect in general. But another approach can be 
used for the study the gauge-fermion systems on the lattice. Indeed, the exponent 
action summed (integrated) over the degrees of freedom of the system is necessary 
only for finding the quantum transition amplitudes. Though in the continuous 
theory the method of finding of transition amplitudes with the help of functional 
integral of the exponent action is very effective, we see that in the lattice theories 



11 



on the such way can arise serious difficulties by reason of lack of the action itself. 
Fundamentally this obstacle can be overcame by means of the direct definition of 
transition amplitudes omitting the stage of action construction |U|. This approach 
seems more perspective when studying the fundamental phenomena, especially in 
the case of discrete space-time: it is not ruled out that exactly quantum transition 
amplitudes are the most elementary mathematical objects in fundamental physics. 
Such approach is realized in this work. 

According to the aforesaid the question of the representation of the partition func- 
tion (quantum vacuum- vacuum transition amplitude) (12.131) in the form of exponent 
action summed over the degrees of freedom of the system is not correct entirely. In- 
deed, the considering sum over closed surfaces is nothing but the partition function 
of the gauge Ising model, that is summed over the system degrees of freedom expo- 
nent lattice action of the gauge field without matter with the gauge group Z^. But 
the sum over surfaces with boundary unlikely can be expressed in a similar manner 
as an integral over fermions of the corresponding exponent, though this sum must 
be interpreted as a Majorana fermions contribution. The last statement follows 
from the comparison of the deciding properties of the weights of the surfaces with 
boundary (the signs of the weights depend on the boundary surface shape only) 
in the considered model and in the case of the lattice Majorana fermions (see Ap- 
pendix B). But the exact coincidence of these weights unlikely can be achieved by 
means of the search of the appropriate lattice fermion model because the fermion 
doubling problem prevents. The correct problem statement seems to be as follows: 
Suppose that in the model presented by the quantum transition amplitudes (but 
not by the action) the second kind phase transition exists, that is the long-wave or 
macroscopic limit takes place. Then one must study the properties of the ampli- 
tudes in the long-wave limit and compare them with the properties of amplitudes 
in the continuous theories . As a result the corresponding identification can be 
determined in the phase transition point. 

Note that the free energy of the studied model has a peculiarity by the parameter 
ijj (see below). This means that there is the phase transition in the model. But the 
nature of this phase transition is not studied here. 



III. THE FIRST STAGE OF DIAGONALIZATION OF THE MATRIX O. 



In the Appendix A the formulas are given which express the partition function 
(12.131) through the eigenvalues of the orthogonal rotation matrix O x>y (see (I A 1 1 ) - 
(I A4[ ) ) . Therefore here and in the subsequent section the problem of finding the 
matrix O x ^ y eigenvalues is solved and the solution is given in (14.171) . The reader 
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wishing to omit these bulky calculations can continue the reading of the paper 
starting from the Eqs. ( 14.171 ). 

In statistical limit, when M, N, L —> oo, the problem of matrix O xy diagonal- 
ization simplifies radically since the matrices O x ,x+z and cj X;X+z become dependent 
only on z at relatively small distances from boundary of the lattice. This property is 
known as translational invariance. Therefore the diagonalization of these matrices is 
performed by means of Fourier transformation, i.e. by passing to quasi-momentum 
representation. The following complete orthonormal set of functions on the lattice 
is used for that purpose: 

\P) = %(rn) = -L e^, \q) = %(n) = -±= e^, \r) = * r (l) = -L e iH , 



A. The first step 



tt(M-2) tt(M-4) 2tt 
P ~ M ' M ' ' ' ' ' ' M' ' ' ' ' 7r ' 

iz{N-2) ir{N-4) 2tt 
Q - ^ , ^ , 0, TT, 

tt(L-2) tt(L-4) 2tt 

r = — 7r 

L ' L ' ""' ' L' '"' ' 

|k) = * k (x) = %(m)^(n)* r (/), k = (p, g, /), 
^*k(x)^(x) = 4k' <— ^ * k (x)* k (xO = (3.1; 

X k 

Let's divide the problem into the successive series of steps. 
It follows from Eqs. (12TT41)-()2TTBT) and (E2D 



Ul h tM,u:-0]i;-X. (3.2) 
= ...O^O^O^ ... . (3.4) 



To find the matrixes O^y^ we use the elementary equations 

(A + £x (A + /i( x (y) = (cos ip% + (sin ^>)£ y , (3.5) 



(A + t4 x £ y y £ y (A + ii£ x Q = (cosip)£ y - (sin^)£ 
(A + ^ X L) & (A + =6, z^x,y 
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(3.6) 
(3.7) 



which follow from (ED) and (ESI). In (l3^ - (l377fl a; ^ y. 

The direct calculation with the help of Eqs. (EEI), (EH) and (EJ3)- (O) gives 
(remember that x = (m, n, /)): 



U \n,l) a ^ U {n,l) = 



(cos 2 ^) a x + (cos 3 ip sin 



x+ei 



+ 



+ (cos 2 ip sin if) 7 x+ei + (- cos ip sin 2 ^) a x+ei +e 3 + (- sin 2 ip) /3 X _ 



-ei+e 3 



+ 



(- cos 2 ip sin 2 ip) a x+ei+ e 2 + (cos ^ sin 2 if) 7 x+ei +e 2 



W; i) = (- sin if) a x - ei + (cos 3 if) /3 X 



+ 



+ 



(- cos 2 ip sin if) a x+e2 + (cos ip sin -0) 7 x+e2 



+ 



(3.8) 



(3.9) 



cos ip sin if) a x _ei + (— cos 2 ip sin 2 ^) /3 X + 



U (n,l)^ U (n,l) = 

+ (cos 3 ?/>) 7 X + (— cos 2 ip sin ip) a x+e3 + (— cos ip sin /? x 



-e 3 



+ 



+ 



(cos ip sin 3 ip) a x+e2 + (- sin 3 ip) 7 X _ 



e 2 



(3.10) 



^(n,0 ax+e 3^( n '0 — 



'sin?/;)7 x + (cos^)a x . 



e 3 



?) /3 x+ e 3 W (nj ;) = (cos ^ sin ^) 7 X + (- sin 2 ip) a x+e3 + (cos ip) 



x+e 3 



(3.11) 



(3.12) 



U ] M a x+e2 U M = (sin^)^ x + (cos^)a x 



-e 2 



(3.13) 



All other ^-matrixes in 



cos ^ sin ip) (3* + (sin 2 ^) a x+e2 + (cos ip) 7 x+e2 . (3.14) 



remain unchangeable. Therefore the right- 
hand sides of Eqs. (13.81 )- ( 13TT41 ) define completely the matrixes O x U y ■ It is evident 
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that if the vector x is not too close to the boundary of the lattice, then Eqs. (13.81) - 
(13.141) are invariant relative to the synchronous change of the number m in x in 
all summands of these equations. In the lattice size infinite limit this means the 
translational invariance relative to the shifts m — ► m + m' along the first axis. The 
translational invariance is used for the partial diagonalization of the matrix O 
by Fourier transformation along the first axis: 



(n,l) 



-P): 



■m 



tnj(p)i (n']l'(p')\ = 2 kj 6 (n,lUri,l') 6 P>P'i ^ 3 = h 2 > 3 , (3.15) 



<,)e!i'(*w)= E 

j,n",l" 



m, n' , I', i; m+m',n", I" ,j 



(3.16) 



The matrix in the square brackets in the right-hand side of Eq. (13.161) is designated 
[0^{p)] n , v vn „ v , . and it is calculated with the help of Eqs. (15^jl-( 1XT4l) . The 
essential part of this matrix is written out below. Under the "essential part" here 
the part of matrix elements with V = I, (/ + 1) and I" = I, (I + 1) is implied. Indeed, 
the linear operator \0^ n,l \p)\ nl y ., n „ v , . acts trivially for the quantities £n,i'{p) with 
V 7^ I, (I + 1). Therefore it is convenient to order the quantities £ n ,i'(p) as follows: 



. . . , v n -i(p), v n (p), v n+1 (p), . . 
where v n (p) are the six-dimensional vectors of the form 



(3.17) 



Vn(p) = 



Pn,l{p) 
Tn,/(P) 

a n ,i+i(p) 

Pn,l+l(p) 
\ln,l+l(p) J 



(3.18) 



In the such designations the matrixes [0^ n,l \p)] n , „ take the form of block-diagonal 
matrixex, and each block is the 6 x 6-matrix, the index n" grows to the right and 
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the index n' grows down. 



M {p) 



\ 



1 

10 

A{p) B{p) 

D CO 

1 



(3.19) 



n'=n+l 



J 



Here A(p), B(p), D, C are 6 x 6-matrixes the same for all matrixes [C?^ n '^(p)] , „. 
The matrixes A(p) and C are placed on main diagonal on the points n and (n + 1). 
correspondingly. According to Eqs. ( 13.81 )- ( jBTT^j ) and ( 13.161 ) 

A(p) = 



(cos 2 -0) (e ip cos 3 ip sin i/A (e ip cos 2 ip sin 0) (— e ip cos ip sin 2 ■0) ( — e zp sin 2 0) 0\ 

-e _ip sin V) (cos 3 0) 

~ ip cos ^ sin i/j) (— cos 2 ip sin 2 0) (cos 3 ip) (— cos 2 -0 sin -0) (— cos ip sin i/>) 

(sinV) (cos^) 

(cos i/> sin 0) (— sin 2 i/>) (cos 0) , 

l) 



B{p) = 



/ {—e w cos 2 sin z 
(— cos 2 i/j sin 
(cos sin 3 



V o 



(e tp cos sin 2 0) \ 

(cos sin 0) 

(- sin 3 ip) 





0/ 



/ (e ip cos sin 2 0) \ 
(cos sin 0) 



sin 






cos 0) 1 



/ 



(3.20) 



(3.21) 



/ 



C = 



COS0) 





( 



V 



sin 






0\ 

1 
(cos0) 
10 
1 
1/ 



(3.22) 
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D = 



/O (sin if;) 0\ 

OOOO 

(-cos?/; sin if;) 





\0 0/ 



(3.23) 



The matrix (13.191) is unitary. This means that 

A(p) [A(p)] t + B(p) [S(p)]t = 1, CCt + D = 1, 
A(p) D* + B(p) C f = 0. 



(3.24) 



The direct check shows that the matrixes (13 . 201) - ( 137231 ) really satisfy the Eqs. (13.241) . 
Thus the first step (finding the matrixes 0^ n ' l \p)) is completed. 



B. The second step 



To realize the second step in which the matrixes (p) are found we must multiply 
the matrixes 0^ n,l \p) at fixed values of the parameters / and p according to the 
rule in the right-hand side of Eq. (13.31 ). Thus 

/ 



/ 



x 



V 



1 
1 


















A(p) B(p) 
D CO 
1 



= < 



n , n 



\ 



\ 



( 



I 



\ 



1 

A{p) B{p) 









D 





1 
1 
1 















A(p) B{p) 
D C 



c 




\ 





1 
1 



X 



7 



(3.25) 



n'=n+l 
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From here we find by means of the direct and simple calculation 

/ 0, for n! < -1 

D, for n' = — 1 



n, n+n 1 



= < 



(3.26) 



CA(p), for n' = 

< C[B(p)} n 'A(p), for ri > 

It is seen from Eq. ( 13.261 ) that the matrix [(9^(p)J , „ depends only on the differ- 
ence (n" — n'), i.e. it is translation invariant (at some distance from the boundary) 
under the shifts along the second axis [12J . This fact enables to perform the further 
partial diagonalization of the matrix , „ by means of Fourier transforma- 

tion along the second axis: 

6(p, q) = *?M£»,i(p) = $(~P> 

n 

£ z (<) (p, q), ^\p', g')] + = 2S itj 6i tV S Ptff S qt ^ i, j = 1, 2, 3. (3.27) 
Now equation (13.31) takes the form 



wj^fa q )u {l) = J2 

j, l" 



£ b (i, (p) 



n, I', i; n+n', I", j 



$ } (p, q)- 



(3.28) 



In (13.281 ) we returned to the accounting of all values of indexes /. Let's introduce 
also the designation 

V \oV\p)} e iqn ' = \o®(p, q)] . (3.29) 

n 

Remind that the matrix (13.291 ) acts trivially for the elements (p, q) with I" ^ 
I, (I + 1). According to ( ED ), 

[5(j9)] n = {->c) n - l B{p), K = sm 2 i/j [e ip cos 2 ^ + sin^] , (3.30) 

Taking into account Eqs. (13.261 ) and (13.301 ) we rewrite the matrix (13.291 ) in the form 



0(0 (p, q ) = \ e -n D + c(l + 



on 



■B(p) A(p) = 



E(p, q) F(p, q) 



H 



G 



(3.3i; 



1 + e^x 

Here E(p, q), F(p, q), H, G are 3x3 matrixes which are calculated with the help 
of Eqs. (Km- (EH}: 

E(p, q) = 



1 + e*«x 



x 
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(3.32) 



F(p, q) = - — — x 
1 + e iq x 



x 



(— e ip cos ip sin 2 ip [e tq sin ip + l] ) 
(e iq cos 2 if) sin 2 ^ [e ip sin ^ - l] ) )((cos^)lo), (3.33) 

sin ^ [cos 2 ip + e ip sin 2 ip{e iq + sin ?/>}] ) 



(3.34) 



(sin^) 

# = | (cos^ sin ^) ] • (3.35) 


The unitarity condition of the matrix ( 13.311 ) is contained in the following equations: 




E(p, q) [E(p, q)]^ + F(p, q) [F(p, g )]t = 1, GG* + H H* = 1, 

E(p, q) H ] + F(p, q) G ] = 0. (3.36) 

The last equalities are verified directly with the help of Eqs. (13 . 321 ) - ( [37351 ) . 

Further it is necessary to return to the accounting of all values of index / and 
at that it is natural to order the quantities &(p, q) as follows (compare with the 
ordering in ( 13.171 ) and ( 13.181 )): 

. . . , wi-i{p, q), wi(p, q), wi +1 (p, q), . . . , (3.37) 

where wi(p : q) are the three-dimensional vectors of the form 

/ ai(p, q) \ 

wi(p, q) = Pi(p, q) ■ (3.38) 
\7z0, q) ) 

In this notations the matrix [(9^(£>, q)\ „, takes the form of a block-diagonal matrix 
with the blocks of the size 3x3, and the index I" grows to the right, the index /' 
grows down, so that the matrixes E(p, q) and G are placed on the main diagonal 
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on the places with indexes / and (/ + 1), correspondingly (compare with ( 13.191 )) 

/ \ 

1 

10 

E(p, q) F(p, q) 



(O (p, q) 



i', I" 



V 



H 




G 
1 



i'=i 
l'=i+i 



(3.39) 



/ 



The second step on finding the matrix 0^(p, q) is completed hereon. 

C. The third step 

On the third, the last step, the matrix O is constructed by means of the ordered 
product of the matrixes O"' according to the rule ( 13.41 ). As a result we obtain the 
matrix invariant relative to the translations along the third axis. The calculations 
are identical to that on the second step, the only difference is in the following 
substitutions: 



A(p) E(p, q), B(p) F(p, q), 

C G, H, 
so that instead of ( 13.261 ) now we have 

0, 



(3.40) 



P(P, q)} 



l,l+l' 



for V < -1 
for V = -1 
for l' = 



(3.41) 



H, 

GE(p, q), 
{G[F(p,q)} l 'E(p,q), for V > 

Using the translational invariance of the matrix ( 13.411 ) we pass to the full Fourier- 
components of ^-matrixes: 

£(p, q,r)=^2 ^r(l)£l(p, q) = £\-p, ~q, -r), 



£®(p, q, rU (i) V, q', r')] = 2<^,A,Ar', i, j = 1, 2, 3. (3.42) 



U^(p, q,r)U = 



£[0(p, q\ l+v e iH ' 



€{p, q, r) 



(3.43) 
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£[0(p, q)] lJ+l ,e irl ' = 0(p, q, 



V 

f . / e ir , \ 1 (e ip + e l<? ) cos 2 ib sin 2 ^ 



1 + e* r cr 



1 + e t( ix 



Here it was used the fact that according to Eq. ( 13.331 ) 

[F(p, grf = (-o^Ffe g). 
3 x 3-matrix (D(p, q, r) in ( 13.441 ) is an unitary matrix: 

0(p, q, r) [<D(p, q, r)] f = 1. 



(3.44) 



(3.45) 



(3.46) 



The last equality is verified easily with the help of the definition ( 13.441 ) and equalities 

(E3SD, (E25D. 

Thus the problem of finding the free energy of the system is reduced to the visible 
problem of finding the eigenvalues of 3 x 3-matrix 0(p, q, r) ( 13.441 ). 

IV. THE CALCULATION OF THE EIGENVALUES OF THE MATRIX 0(p, q, r) 



Let's give the products of the matrixes in ( 13.441 ) in the evident form: 
GE(p, q) = - — — — x 



1 + e^x 

(cos 4 tp) (cos ip sin tp[e zp cos 2 ip + e~ iq + sin 3 ip\\ 

X I (y~ cos V' sin V 1 |^cos 2 ip sin ip + e Lq cos 2 ip -J- e*^ — ^ 4 " 9 ^ sin i/> + e _ip j (cos 4 -0 — sin 3 i/![e ip cos 2 i/i + e -19 -J- sin 3 0]) 
(e~ ip cos 2 ip sin 0[e* p sin -0 — 1]) (— cos ip sin 0[e — %q + sin -0]) 

(e ip cos 4 ip sin i/> [e 19 sin ijj + 1]) 



(-cos 3 ip sin [e i(p+g) sin ip + e ip sin 2 V - e iq cos 2 V']) J ; (4.1) 
(cos 2 ip [cos 2 i/; + e zp sin 2 ^{e* 9 + sin ) 



GF(p, «)£(p, q) = 



(1 + e i9 x) 



e ip cos 2 V> sin 2 V \e iq sin ^ + ll 



\ 



X(Ul, U2, U3j, 



"1 



(cos ip sin 2 V [e l(p+9) sin ip + e ip sin 2 t/j - e iq cos 2 v] ) 
y (— sin ip [cos 2 ip + e tp sin 2 ip{e lq + sin?/)}] ) y 

li2, W3) = ( (cos 4 ip — sin ^ |V 9 cos 2 ip + e l ^ p+IJ ' sin tp + e~ lp J ) , (cos ip [cos 4 ?/> + sin ip{e lp cos 2 ^ + e~ %q + sin ?/;}] ) , 

([ e <p + e <9] cos 4 V sin ^) J. (4.2) 
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Note that all columns of the matrix (14.21) are proportional to the third column of 
the matrix (14.11 ). 
For completeness let's give also the following formulae: 

0(p, q, r) = \e~ ir H + G (l + Y^~ F ^ <?)) E ^ «)} > 
(e ip + e zq ) cos 2 ip sin 2 ip 



a = 



1 + e^x 

X = [l + e tq K)[l + e ir a) = 1 + e iq sin 3 i/j + 



x = sin 2 ip [e ip cos 2 ^ + sin ip] 



e i(p+q) e i(p+r) ^(g+r) 



/0 (sin0) 
if = I (cos?/; simp) 
\0 

We must solve the cubic equation 

det ||C> - p^H = 
relative to the variable p. This equation is rewritten as 

p 3 - (trO)p 2 + (mn + rn 22 + rn 33 ) p - det O = 0. 
Here is the minor of the matrix element 0{j. 



cos f/>sin -0, 
(4.3) 

(4.4) 



(4.5) 



(4.6) 



With the help of Eqs. (14TTT)-(1Q) we find: 



(tr O) = 1 - -, 77 = (1 — 3 cos 4 iP + sin 6 0) + sin 3 (e^ + e^ 9 ) = fj. (4.7) 



Since the coefficients in Eq. (14.61 ) are the invariants which do not depend on the 
matrix O representation, we have: 

(mn + m 2 2 + m 33 ) = (pip 2 + pips + P2P3) = 

= (Pi 1 + P2 1 + Pz 1 ) P1P2P3 = (tr C f ) det O = (tFO) det O. (4.8) 

Here it was taken into account that O^O = 1. Thus it remains to calculate the 
quantity det O only. 
Let V4, i = 1, 2, 3 be the vector-columns, so that 

1 



1 + e^x 



(Ui, V.2, V.3) . 



(4.9) 
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Then the matrix O is represented as follows 



O 



v.i 



V.2 



, W.3 



e ir sin ?/> 



cos 2 ip (1 + e i9 x) x 



U.3 X (ui, U2, 0) 



1 



w.3 = — -u.3 + r, r = e ir sin?/; cos?/; 



(4.10) 







and (wi, W2, u%) is the row matrix in (14.21) . 

Let's add to the first two columns of the matrix O such columns, proportional 
to u>3, which cancel the second summand in Eq. (14.101 ). Thus the matrix O' is 
obtained with the same determinant: 

e ir sin ip 



AO = 



a = o + ao 



sin 2 ip 



cos 2 1\) (1 + e^x) 

/ V.i V. 2 



W.3 X (u\, U 2: 0) , 

, w. 3 ) + 



1 + e^x 1 + e^x 



1 



cos 2 "0(1 + e ig x) 



cos?/; ] x w 2 , 0) 




(4.n; 



Now let's subtract the first row of the matrix O' multiplied by (cos?/;) from the 
second row of this matrix. As a result we obtain the matrix O" the determinant of 
which coincides with the determinant of matrix O: 



o" = 



/ 
\ 



JH 4 Tjj-eiq cos 2 j, 5in 3 ^,_ e '(-P+9) sin 4 ^_ e -'P sin 3 ^ 
cos 2 ,/.(l + e'9«) 

^(™ ! </>+e* 9 cos 2 i> sin^ + e't-P + s) sin 2 V^e - *" sin 
cos^ ^isini/j^sini/i- e~^ p ^ 



:in*(=in^ + e i P cos 2 * + e" i<! ) 
co S +(l + t>5x) 



lp cos 2 ^ sin j>- 



i^sinV.(sin^ + c~ i ' J ) 



cos 4 if S i n +(t i <P+«) sin^ + c i P)+c- ir X sinY' 
1 i, sin^(c i (P + l) sin^(l+cos 2 ^J + e'P-c'? cos 



(4.12) 



: t( C DB 2 ^ + e i (P+9) sin 2 ^. + e*P sin 3 </>) / 
X / 



det G" = tetG. 
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Let m'-A be the minor of the matrix element O"--. We need the following minors: 

cos 2 ip sin ip {eT %v + e~* 9 ) 



m 13 = 



1 + e l iH 

.. e~ ip sin 3 i/j — e~ iq cos 2 i/j sin if; + e - *^"^ sin 2 ^ 
23 cos^(l + e^x) ' 

„ cos 2 ip + e~ ip sin 3 ^ + e ~ l(:p+9) sin 2 tp 
77133 ~ cos 2 ^(l + e^x) ' ^ ' ' 

To compute the determinant let us decompose det O" in the last column of the 
matrix O" . Thus we obtain: 

detC = ^. (4.14) 

X 

With the help of formulae ( 14.71 ) and ( 14.81 ) we obtain: 

(ran + m 22 + m 33 ) = - — -• (4.15) 

X 

As a result the equation ( 14.61 ) takes the form: 

f-(l-*)? + h-l) P -* = (p-l)fo + lp + *)=0. (4.16) 
V XJ \X XJ X V X XJ 

Now we find easily all eigenvalues of the matrix O: 



Pi,2\p,q,r) = - , P3 = l- 4.17 

Note that the expression (4%X ~~ 7 ? 2 ) which is under the sign of square root in ( 14.171 ) 
is a non-negative number. Only under this condition we have 

|pi,2,s(p, q, r)\ = 1. 

The last equalities follow also from the unitarity condition of the matrix 0(p, q, r). 
Otherwise it would be 

|pi,2(p, 9, r)\ ^ 1. 

The last inequalities mean the violation of the matrix q, r) unitarity. 
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V. THE FREE ENERGY 



According to (Elj) and ( 14171 ) 

tr U = <2MNL/2 [] + ^) + v'S) = 2™ /2 I] * + * " ^ 



a/ yv 

= 2 AfiVi/2 TT (1 + 3 cos 4 ^ - sin 6 ^) + 2i/(p, g, r) cos 2 ^ sin 2 ^ 

p,q,r v 

z/(p, g, r) = [cos(j9 + g) + cos(]9 + r) + cos(g + r)} . (5.1) 

The free energy is of interest. Up to inessential summand the free energy has the 
form 

MNL 



T = -T\nZ = T- 



167T 3 



/TT pTT HIT 

dp dq dr(lnx + lnx) 
■TT J —TT J —TT 



MNL f n f n f 77 ( 1 

-T - 3 / dp dq drln I (l + 3 cos 4 ?/; - sin 6 ?/;) + 2z/ cos 2 ^ sin 2 ^ >. 

(5.2) 

Here it vas taken into account that at M, iV, L — > oo the substitution 

p q r J-TT 17 _7r 

is valid. 

The first summand in the right-hand side of Eq. ( 15.21 ) is equal to zero. This fact 
is the consequence of the slack inequalities 

W\ < 1, M < 1, (5.3) 

so that the equalities in ( 15.31 ) take place only on the subset of the zeroth measure 
in the space of the variables {p, q, r}. Indeed, according to ( 14.31) and ( 15.31) we have 
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the following chain of equalities: 



»7T />7T 



dp / dg / drlnx = 



t/ — 7T t/ — 7T 

''IT PIT PTT ( 



dp dq dr ln(l + e !,| x)+ln(l + e l V) 



— 7T 1/ — 7T 1/ — 7T 



•>7T /»7T 



= E L ^— d(? / df|e in W,p)r + e-[* M )r| = 0. 



n=l 



-7T <V — 7T 



(5.4) 

Prove the slack inequalities ( 15.31 ). 

According to ( 13.451 ) F 2 = —aF. Let's take the modulo of the matrix element 
1, 1) of the last equality: 



kll*i] 



< 



\ 



= \/ (F F^)n (F^F)iv (5-5) 
We have also from the unitarity condition ( 13.361 ) and the evident form of the matrix 

G QSM) 

(F f F) ! ! = 1 - (G^G) 1 1 = cos 2 ip sin 2 ^. 
By means of the direct calculation we find also that 



(5.6) 



1 + 2 
11 cos z ip 



(5.7) 



Combining all formulae (15. 51 ) -( 15771 ) we obtain: 

\a\ < J sin 2 ip(l + cos 2 ip) < 1 
Using the evident expression ( 14.31 ) for x we find: 

\x\ < cos 2 ip sin 2 ip + sin 3 ip < 1 for < < 7r. 
Thus, both inequalities ( 15.31 ) are verified, and therefore 



T=-T—^\ dp dq drx 



8tt 3 



-7T <V — 7T 



X 



In | (1 - sin 2 ^) [sin 4 ^ + 2(i/ - 1) sin 2 ip + 4] 



(5.8) 
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It is seen from the expression for free energy (15.81) that the earlier separation of 
the second axis in the previous formulae for the intermediate quantities, including 
the eigenvalues of the matrix 0(p, q, r), was an artifact: the free energy of the 
system, which presents the physical interest, is completely symmetric relative to 
mutual substitutions of the all three axes. 

From the definition of the parameter ip (see the second point) it is evident that 
this parameter plays a part of the temperature T. At T — > oo the angle ip — > 0. 
This follows from the fact that at T — ► oo the contribution of surfaces into the 
partition function ( 12.131) tends to zero. That corresponds to the high-temperature 
limit in the gauge compact lattice theories. It is also seen from the comparison of 
the partition function (12.131) with the corresponding quantity of the lattice gauge 
theory (with the gauge group Z<i) that the decrease of the temperature from the 
infinity up to zero means the surgeless increase of the angular parameter ip from 
zero up to (tt/2): 



The answer for the important question is need: do the phase transitions exist at 
the temperature reduction from the infinity up to zero, and what is its kind? 

Since the phase transition point coincides with the peculiarity of a free energy by 
the temperature, so one needs look for that values of ip at which the expression 
(15.81) can have a peculiarity. Just as in the case of the two-dimensional Ising model, 
the only possibility to have a peculiarity in the free energy consists in nullification 
of argument of the logarithm situated under the sign of the integral over quasi- 
momentum. This can occur at some values of the temperature and quasi-momenta. 
Then the integral over quasi-momenta near the peculiarity gives irregular contribu- 
tion into free tnergy and other thermodynamic quantities. In the case of the two- 
dimensional Ising model there is the only point in the three-dimensional manifold 
of the totality variables "two quasi-momenta+temperature" in which the logarithm 
argument vanish: p = q = and T = T c . 

In the considered three-dimensional theory the situation differs qualitatively from 
the indicated one in two-dimensions: for ip = (T = oo) the logarithm argument 
in (15.81) is positive always; for p = q = r = the logarithm argument in (15.81 ) 
vanish only for ip = 7r/2 (T = 0). But there is a possibility of the logarithm 
argument nullification in (15.81 ) at finite temperatures < ip < tt/2 and nonzero 
quasi-momenta. Indeed, at the quasi-momentum variations the parameter v varies 
in the following range: 



dT 



<0, 




(5.9) 



3 < v < 3. 



(5.10) 
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To zero the logarithm argument in (15.81 ). it is necessary realization of the condition 

x\ + %y - l)x c + 4 = 0, x = s\u 2 %Ij. (5.11) 

Only the real solutions on the segment < x c < 1 are of interest. This is possible 
only under the condition 

-3<z/<-5 (5.12) 

The only solution has the form 

x c (v) = (1 - v) - v / (l-^) 2 -4, (5.13) 



at that \j x c {y) is a monotonic increasing function on the segment (15.121) . so that 



y/x c (-S) = sin^ c = ^4- Vl2 < 1 (5.14) 
is its minimum value, and 

v/x c (-3/2) = sinV c = l (5.15) 

is the maximum possible value. 

The point v = —3 is the evolved point in the space of quasi-momenta since it is 
obtained for two values of quasi-momenta p, q, r only. Indeed, for that the quasi- 
momenta must satisfy at least one of the following eight equations (all signs "plus" 
and "minus" in the right hand sides of Eqs. ( 15.161 ) are mutually independent): 

p + q = ±7r, p + r = ±7r, q + r = zbrr. (5.16) 

Since the system of linear equations ( 15.161 ) is non-degenerate, so for each value of 
the right hand side in ( 15.161 ) there is only one solution for quasi-momenta p, q, r. 
Taking into account the fact that the quasi-momenta are defined modulo 2rr we can 
reduce all eight solutions of the systems ( 15.161 ) to two independent solutions: 

M) = M) = r (i) = 1 

fc 1 c 2 ' 

pf) = q f) = r P) = -I. (5.17) 

On the contrary, for all values of v from the half-interval 

-3 < v(p, q, r) < ~ (5.18) 

there is a whole continuum of the quasi-momenta values for which the logarithm 
argument in ( 15.81 ) vanish. 

The further study of the free energy properties and the properties of a phase 
transition in the considered model will be performed in an another work. 
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VI. CONCLUSION 



In the present work the calculation method of the partition function of the 
gauge system with the gauge group Z<i coupled with Majorana field on the three- 
dimensional cubic lattice is suggested. Actually the sum over closed surfaces and 
the surfaces with boundary is computed. The surfaces can be with self-itersections 
and their weights are proportional to the factors fi 2S , where S is the number of 
the faces of the surface. The weights of the closed surfaces are positive always. 
The total sign of the surface with boundary depends essentially on the boundary 
configuration but not on the surface form with the given boundary. This sum is 
computed completely, it is represented as a threefold integral over quasi-momenta 
and it is the function of one parameter (the coupling constant or the temperature). 

It should be noted that there is little examples of the integrable and "rational" 
(that is having physical interpretation) statistical systems in three dimensions, show- 
ing a phase transition relative to the temperature. Thereupon I want pay attention 
to the work of A.B. Zamolodchikow [7f| in which the Yang-Baxter triangle equa- 
tion has been generalized to the three dimensional case where the corresponding 
equation is named as tetrahedron equation. 



APPENDIX A 



It is known [3j that the rotation matrix in spinor representation (12.141) can be 
expressed in the form 



U = ex P Q^x,y7x7yJ , ^x,y = -W y , x , 



and 



U^JA = C x , y 7 y , C> x , y = (e w ) X)y = £ x , y + u Xiy + i uj x . 



(Al) 



(A2) 



2j ~x,z^z,y 

The trace of the matrix IA is expressed simply through the eigenvalues of real or- 
thogonal matrix O x ,y Let the set of numbers 

(Pi, Pi, f>2, fo, ■ ■ ■ , P3MNL/2, P3MNL/2) = {Pk, Pk}, k = 1, . . . , 3MNL/2 (A3) 
form the complete set of eigenvalues of the matrix x . y . Then (see Appendix A) 

3MNL/2 r , . -. 3MNL/2 f / , \ n 3MNL/2 



tiU 



n 

k=i 



2ch 



In Pk 



n 

k=l 



Pk 



2 cos 



(f>k 



n (v^+v^)' 

(A4) 
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Let us prove the formula ( 1A4I ) . 

Let lv x k \ v x k ^\ , k = 1, . . . , 3MNL/2, be the complete orthonormal set of eigen- 
vectors of the matrix X}V , so that the eigenvalue p k (p k ) corresponds to the eigen- 
vector y^J • Further also the designation 

{v x k \ ^ fc) } = K}, o = l, ...,3MNL 

is used. We shall consider the introduced vectors as vector-columns and the upper 
indices T and ' denote the transposition and Hermitian conjugation of vectors and 
matrices. By definition 

v (k)T v (k') ^^vi k H k '\ vWvW^Y^vM. (A5) 

X X 

The given definitions imply the following formulas: 

««V = 0, v^vW = 5 kk > : (A6) 

S xa = v a x oiS = (y®, v( 2 \ . . .) , = 5 abl (A7) 

(tfOS)^ = diag (pi, p 1? p 2 , p 2 , . . . ) = D ab . (A8) 
It is shown in [3] that 

(SfluS) ab = diag (In pi, -In pi, lnp 2 , -lnp 2 , . . .) = A ab . (A9) 
Due to (Ell) and (011) we have 

E \ ^ a = E E 3 fe 5 -) A ^ (^y&) • ( A1 °) 

£,2/ x,y a,b 



2 3MNL/2 x 2 3 ™/ 2 -matrixes 



4 = 7^, c fc = 7x^ } (All) 

possess all properties of fermion creation and annihilation operators. Indeed, in 
consequence of (12.11 ) and (IA6f ) 

[°k: 4']+ = <w> [c/t, c k >] + = [4, 4']+ = °- ( Al2 ) 

According to the definitions (IA7[ ) and (I A lip we have 



(A13) 
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With the help of Eqs. (E2), ( 1X121 ) and ( 1XT3I ) the quantity ( IXTUj ) is rewritten as 

3MNL/2 3MNL/2 



x,y k=l k=l 



(In p k )c\c k - ^lnp k 



1 



(A14) 

Equality (IA4[ ) follows immediately from (IA14H since the calculation of the trace 
in terms of ^-matrixes is equivalent to the calculation of trace in terms of the 
corresponding fermionic operators ( lAllf ). 



APPENDIX B 



Let us consider he Majorana spinors on the three-dimensional cubic lattice an 
their contribution into the partition function. 

In the simplest case the action of the Dirac fermions on the cubic lattice has the 
form 

3 



7 1 = ^, 7 2 = ^ y , 7 3 = ^, Ul ei U x , ei = l. (Bl) 

Here Ui are the Pauli matrixes and U Xjei is the gauge field. The Fermi-fields and 
are the elements of the Grassman algebra and all their elements are considered as 
a mutually independent variables. The Fermi contribution to the partition function 
is defined as the integral 



z d{U} = Y[ J d*p x dip x ex.pS Dl 



(B2) 



where d ^x an d d ^x denote the products of the differentials of the both components 
of the corresponding spinors. 
The Majorana spinors are determined by the following system of identifications: 

^x = ~€l 2 - (B3) 
The "electrical current" of the Majorana spinors is equal to zero identically: 



= Oxl, ^2) 



-i J H 1 J ' [ i j ' v -1 



^xl 
^x2 



- "Wxl + Vi), "^xlfe + fe^xl)} {0, 0, 0} 



(B4) 
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The last identity in ( 1B4[ ) follows from the fact that all components of the field ip are 



odd elements of the Grassman algebra. The majorana action is obtained from the 
Dirac action (IBID by means of the substitutions — ► — an d division by 2: 

3 

S M = ~\ E E (Vx,e^x+e t " ^-e, e^x-e.) • (B5) 

x i=l 

The contribution of Majorana spinors to the partition function is determined by the 
Grassman integral 

Zm{U] = II / d^x 2 d^xiexp<S M . (B6) 

X ^ 

Since 

/ d ^x2 d Ipxl (^xlfe)™ = 



1, for n = 1 
0, for n / 1 



so 

d^x = 0, 



/ 

-« y d^x- (V'xV'x) = ~ i J d^x2d^xl 



^xl^A _ 2 
fe^xl J ~ 7 ' 



/ d^x-(^x)®(^)=0, (B7) 

and so on. 

Draw on the lattice the system of closed broken contours without intersections and 
self-intersections. The elementary link of each contour is an edge [ X) j (see Section 
2.2) and each vertex belong to one and only one contour. Each contour is oriented, 
that is the direction (arrow) is assigned to each edge of the contour, so that the 
continuous movement along the arrows reduce to the whole round of the contour. 
We shall call the edge as positive (negative) oriented if its arrow is directed in the 
line of the positive (negative) direction of the corresponding lattice axis. Below the 
particular case of the elementary closed contour based on only one edge with both 
orientations at once is described. 

The nonzero contribution in the integral (IB6I) give only those summands in the 
exponent expansion in the quantities 



^7V£Wx+e. t ) (B8) 
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and 



(B9) 



in which these quantities are in the first power. The factor (IB8I ) corresponds to 
each positively oriented edge and the factor ( 1B9I ) corresponds to each negatively 
oriented edge [131 ] . The both factors ( IB 8ft and (IB9I ) correspond to the elementary 
closed contour based on the edge [ x ^. Thereby, under the sign of the integral in each 
vertex x there is the element of the Grassman algebra (Vx^x) in the first power 
giving a nonzero factor according to (IB 7ft . Note that the linkage of this rule with 
the contours orientation in the Majorana case is due to the interaction between 
fermions and gauge field. Indeed, the calculation of the integral over the gauge 
field in the frame of the high temperature expansion leads to the nullification of all 
contributions from the non-oriented closed contours (in the case of the higher gauge 
symmetry than Z^) which do not vanish as a result of the fermion integration (IB 6ft . 
Here the specific question is interesting for us. Therefore further we put t/ x ,e l = 1- 
From the aforesaid and with the help of Eqs. ( 1B5I ). ( 1B7I ). ( 1B8I ) and ( 1B9I ) the 
following rules follow for the calculation of the integral (IB6I ): 

1) Let's draw on the lattice the system of the outlined above closed broken oriented 
lines and, rounding each contour successively along the arrows, relate successively 
to each edge [ X) j the factor (1/4) y in the case of the positive orientation of the edge 
and the factor (—1/4)7* m the case of the negative orientation of the edge. 

2) After ending the rounding process of each contour let's calculate the trace of 
the ordered product of the 7-matrixes corresponding to the contour according to 
the rule 1, and add the factor (—1), common for the whole contour. The obtained 
number is called as the factor of the contour. 

3) Let's multiply the factors of all contours. The obtained number is called as the 
factor of the system of contours. 

4) For finding the integral (IB6I ) it is necessary to summarize the factors of all 
possible systems of contour. 

Now let's show the change of the sign of the contour system factors on the simplest 
examples. Since only the sign is interesting for us, the others positive factors are 
ignored. 

Let's consider the factor of the elementary closed contour based on the edge [ x ^. 
This contour is represented in the Fig. 4. 
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FIG. 4 



According to the given rules the factor of this contour is 



$ x ,i = (-i)try(-y) = i. 



;bio) 



One of the possible and at the same time simplest system of contours is represented 
in the Fig. 5. In Fig. 5 one of the mutually parallel planes containing the base 
vector ej is represented. Thus, all these planes contain the identical configurations 
of elementary closed contours and each vertex of the lattice belong to one and only 
one contour. It is evident that due to (IB 101 ) the factor corresponding to this system 
of contours is equal to the product of units and thus it is equal to the unity. 
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FIG. 5 



It is supposed that in the subsequent examples of the closed contour systems 
almost all closed contours are elementary, so that the change of the number of 
the elementary contours does not affect on the total sign of the contour system. 
Therefore only that closed contours will be considered and represented in the figures 
which can effect on the total sign of the contour system factor. 

Let's consider the factor corresponding to the closed contour in the plane (1, 2) 
which is represented in the Fig. 6. 
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x + e 2 x + ei + e 2 




x + ei 



FIG. 6 



The factor corresponding to the contour in the Fig. 6 is 



l)tr 7 V(-7 1 )(-7^) = l- 



x,ei,e 2 



(Bll) 



Therefore the total sign of the contour system factor corresponding to the Fig. 6 is 
also equal to unity. 

Now let's consider the factor corresponding to the closed contour in the plane 
(1, 2) and represented in Fig. 7. 

x + e 2 x + ei + e^x + 2e x + 2e 2 



x + e-2 



■x + 2ei + e 2 



x + ei x + 2ei 



FIG. 7 



,2ei,2e 2 



= (— 1) tr7 1 7 1 7 2 7 2 (— 7 1 )(— 7 1 )(— 7 2 )(— 7 2 ) = -1. 



(B12) 



From here it is seen that the total sign of the contour system factor corresponding 
to the Fig. 7 is negative. 

Thus we see that the contour system factor, the sum of which defines the integral 
(IB6I ) . can have either positive or negative sign. It is important that the change 
of the contour configurations leads, generally speaking, to the change of the factor 
sign. 

In conclusion note that the fermion Dirac contribution into the partition function 
(IB2f ) is expressed also as a sum of the contour system factor. But in the Dirac case 
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the sign of the factor does not depend on the contour configurations. It is easy to 
see that in the long-wave continuous limit in a weak gauge field. In this limit the 
lattice action (IB II) transforms into the usual Dirac action 



S D = J d 3 xif> (ij'di - eYAi) i/>. (B13) 

The corresponding contribution into the partition function can be writthen in the 
form 

det - e^Ai) = Const • exp {trln [l - e^d^y'Aj] } = 

= Const • exp \^-e 2 J d 3 xd 3 yti^d^-^A^y^dk^^A^x) - . . .J . 

(B14) 

In the last expansion in the exponent under integrals the enough smooth function 
(H^jOxu i s present. Since the space integrations mean the contour variations, so 
it is seen that the contour variations is not conjugated with the sign variation of 
the corresponding factor. The variations of the contour factor signs of the lattice 
Majorana fermions mean that in the continuous long- wave limit the contributions 
into the partition function from these contours are cancelled mutually. From here 
the impossibility of Majorana fermions- Abelian gauge field interaction in the case 
of continuous theory is seen. This state also follows directly from (IB4I ) . 
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